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SUMMARY 

The construction of interferograms, schlieren, and shadowgraphs from computed flowfield solu- 
tions permits one-to-one comparisons of computed and experimental results. A method for construct- 
ing these images from both ideal- and real-gas, two- and three-dimensional computed flowfields is 
described. The computational grids can be structured or unstructured, and multiple grids are an 
option. Constructed images are shown for several types of computed flows including nozzle, wake, 
and reacting flows; comparisons to experimental images are also shown. In addition, the sensitivity 
of these images to errors in the flowfield solution is demonstrated, and the constructed images can be 
used to identify problem areas in the computations. 




INTRODUCTION 

The development of techniques for comparing computed flowfield solutions with experimental 
interferograms, schlieren, and shadowgraphs is vital for validation of ideal- and real-gas computa- 
tional fluid dynamics (CFD) codes. Infinite- and finite-fringe interferograms can be transformed into 
density fields only for two-dimensional (2-D) and axially symmetric non-reacting flows. In experi- 
mental schlieren and shadowgraphs, several flow features including shocks, shear layers, and 
expansion fans are recorded. When flow solutions are compared to these images, contour plots for 
only one plane of data, usually the symmetry plane, are typically used. If the flow is three dimen- 
sional or the model is free to roll, no single computational plane provides all the information 
necessary for realistic comparisons to experimental images. 

The information required for the construction of interferograms, schlieren, and shadowgraphs is 
contained in the flowfield solutions. By developing routines to evaluate and integrate appropriate 
functions of the refractive index (n), direct comparisons of computed and experimental results can be 
made. Constructed interferograms provide the basis for quantitative comparisons of line integrals of 
n for reacting flows and the integrated density for non-reacting flows. With constructed schlieren and 
shadowgraphs, the computed and experimental locations of shocks, shear layers, expansion fans, 
separation, and reattachment can be compared. Constructed schlieren can also be used to compare 
signs of line integrals of the gradient of n. 

In this paper, CISS (Constructed Interferograms, Schlieren, and Shadowgraphs) is described. 
CISS is software that constructs images from ideal- and real-gas, 2-D, axially symmetric, and 
3-D computed flowfield solutions. The computational grids can be structured or unstructured, and 
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multiple grids are an option. Results are shown for several types of flowfields, and the effects of grid 
resolution and solution convergence are discussed. Comparisons to experimental images are also 
included. 


CONSTRUCTING NUMERICAL INTERFEROGRAMS, SCHLIEREN, 

AND SHADOWGRAPHS 

Construction of interferograms, schlieren, and shadowgraphs from flowfield solutions requires 
three steps: (1) identifying and evaluating the appropriate functions of n, (2) integrating these func- 
tions along lines of sight, and (3) post-processing the integrals to give the desired images. In this 
section, these three processes will be described. 


Functions of the Refractive Index used in Constructing Images 

As light passes through the flowfield, the phase shift and angle of deflection are obtained by 
integrating functions of the refractive index along lines of sight. For ideal and non-reacting gases, the 
refractive index is simply 


n = 1 + Kp 

where k is the appropriate Gladstone-Dale constant for the gas and p is the density. For both equi- 
librium and nonequilibrium real-gas computations, the species mass fractions vary throughout the 
flow, and the refractive index is given by the sum of the contribution from each gas (ref. 1), 

n = l + £ KjPi 
i 

Here, K, is the Gladstone-Dale constant and p, is the partial density for the ith species. 

Interferograms measure the phase shift between reference and object beams caused by variations 
in the refractive index throughout the flowfield. When calculating interferograms from computa- 
tional flowfields, the phase shift is obtained by integrating 

f(n) = ^(n-n 0 ) 

along a line of sight. Here, no is the refractive index of the undisturbed flow and X is the wavelength 
of the light. 

Schlieren are created in the laboratory by passing collimated light through the flowfield, focusing 
this light to a finite-sized point, using a knife edge to block out a portion of the light, and then 
re-focusing the remaining light onto an image plane (fig. 1). The amount of light blocked by the 
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Figure 1. Schematic for constructing schlieren. 


knife edge is proportional to the angular deflection of the light; the deflection is caused by gradients 
in the refractive index. The incremental change in the angular deflection at any point in the flowfield 
is 
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In shadowgraphs, the image plane is not focused in the flowfield, but at some distance from it, 
and the dark and light regions are caused by the concentration and divergence of light. In many text 
books (e.g., refs. 2 and 3), the function used for constructing shadowgraphs is the line integral of the 
second derivative of the refractive index. In CISS, the function used for constructing shadowgraphs 
has been chosen to be equal to that used for schlieren. At the exit of the computational flowfield, the 
angular deflection is known. If the distance to the image plane is also known, then the location on 
the image plane where the deflected light falls is known. The shadowgraph is constructed by adding 
the contribution of the deflected light beams at each point on the image plane. Using this method, the 
thickness of the dark and bright regions is partially controlled by the position of the image plane. 
This dependency of the dark and bright regions on the image plane’s position is observed in 
experimental shadowgraphs. 


Integration Scheme 

Tracing the actual light path as it bends through the computational flowfield and integrating the 
appropriate function of the refractive index along this path is computationally expensive. The 
resources required for producing constructed images can be reduced by approximating the light’s 
path by a straight line perpendicular to the image plane. For regions of the flow without shocks, this 
approximation should have a minimal effect on the constructed images. In shock regions, the index 
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of refraction changes rapidly, and the approximation may introduce errors; however, these errors are 
in many cases on the order of or smaller than the solution and experimental errors. 

This approximation minimizes the need for search routines, and interferograms, schlieren, and 
shadowgraphs can be constructed without any ray tracing. For instance, consider an image defined by 
an M x N array of light ‘cells.’ The value of the integral for the m,n element of the array is given by 

z 

A m,n= J f [ n ( x m>y n >0] d C 
0 

where the image plane is located at z and (x m , y n ) are the x,y-coordinates of the light cell at the 
m n location in the array. The integration path can be broken into several segments 

N z i+l 

A m,n = ^^ J f [ n ( x m’yn»0] d C 

i=l Zj 

In this expression, the line of integration enters the ith computational cell at z\ and exits it at z\+\ 

(fig. 2). In the actual solution procedure, no ray tracing of the line of integration through the flow- 
field is necessary. By assuming minimal deflection of the light through the flowfield and approxi- 
mating the light’s path by a straight line, the contribution of each computational cell to the line 
integral can be found independently, and the order of evaluation and summation has no effect on the 
final result. Therefore, for each computational cell, the integrals 

z i+l 

J f[n(x m ,y n ,0] d C 

z i 

are evaluated for every x m , yn that fall within the cell, and the value of the integral is added to the 
appropriate location on the image plane. 

The test for determining if the point x m , y n falls within the computational cell is as follows. First, 
to minimize the number of tests required, the maximum and minimum x and y values for the compu- 
tational cell are found. Then the surfaces of the computational cell are described by triangles, and for 
each triangle, three cross products are formed: 


(xi-x mn )x(xi-x 2 ) 


(x 2 -x m n)x(x 2 -x 3 ) 


(X3-Xmn)x(x 3 -Xi) 
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Figure 2. Schematic for integration process. 

In these equations, x mn gives the position of the light cell, and xj, X 2 , and x 3 are the projections 
onto the image plane of the vertices of the triangle. Testing is performed only for the x m and y n that 
lie within the minimum and maximum x and y values for the computational cell. The sign for all 
three cross products is the same if and only if the point lies within the triangle, and the triangular sur- 
faces through which the line of integration enters and exits the computational cell can be determined. 
Once these surfaces are determined, interpolated values for z\, z\+\, f(x m , y n , z;), and f(x m , y n , zj+i) at 
the entrance and exit are obtained, and the line integral is evaluated. In CISS, three functions of the 
refractive index are integrated simultaneously: n— 1, (l/n)dn/8x, and (l/n)8n/3y. 

Construction of the Images 

For infinite-fringe interferograms, the image intensity is simply proportional to the sine of the 
phase shift <|>: 

z 

^ = xj (n-n 0 )dC 
0 

For finite-fringe interferograms, an additional phase is added: 

5<|) = k x x + k y y 

and the spacing and orientation of the freestream fringes are determined by k x and k y . 
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For schlieren, the intensity is a function of the angular deflection. As an example, consider a 
horizontal knife edge (fig. 1) and a rectangular section of the initial light source with height h. As the 
rectangular section of the light source passes through the flowfield, it is deflected and exits the flow 
at an angle. When the light exiting from the flowfield is then focussed to a finite-sized dot at the 
plane of the knif e edge, the rectangular element has a height of h', and it has been deflected vertically 
by KE y ; k is a function of the distance between the flowfield and the knife edge. There is also a 
horizontal deflection; however, it is unimportant in this case. If the horizontal knife edge is placed at 
the center of the plane, the ratio of light not blocked by the knife edge to the amount of the original 
light is 


0 

for 

K£ y< 1 
h' s 2 

1 

for 

Ke y > 1 
h’ ~2 

1 K£ x 

2 + h’ 

for 

1 K£ y i 

-i < ir < i 


The intensity patterns in the schlieren are proportional to the square of this ratio. 

As mentioned previously, shadowgraphs are not constructed in CISS from line integrals of sec- 
ond derivatives of the refractive index. Instead they are given by variations in the intensity caused by 
the deflection of the light. The location and size of a square element of light after it passes through 
the flowfield is determined by the deflection of light sources that define the comers of the element. 
The total angular deflection for these light sources is a function of the angle at which the light exits 
the flowfield; that is, shadowgraphs are also governed by e x and e y . The deflected square element is 
then projected onto the image plane, and the amplitude of the electric field at each i,j -element in the 
image plane is increased by 


B f nA ij 
" E o B; 


Here, E 0 is the amplitude of the initial field, Bf n Ay is the area of the intersection of the deflected 
element with the i,j-element of the image plane, and B; is the initial area of the deflected element. 
The intensity pattern observed in the shadowgraph is proportional to the square of the Ey's. 

The time requirements for constructing these images from simulated flowfields are linearly pro- 
portional to the grid size and the number of pixels defining the image. For 2-D flowfields, the image 
construction takes only a few seconds on a CRAY-YMP. For 3-D images defined by 400,000 pixels, 
construction using flowfield solutions with 250,000 grid points takes a few minutes. 
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CONSTRUCTED IMAGES 


Representative images are shown in this section for several types of hypersonic flows including 
ramp, nozzle, reacting, and blunt-body flows. In figure 3, constructed and experimental schlieren are 
shown for a 2-D, 15° compression comer; a horizontal knife edge was used for both the schlieren. 
The shocks, shear layers, and separation region observed in the experimental schlieren are simulated 
by the constructed schlieren. The dark regions in the constructed schlieren are caused by rapid 
increases in density across the computed shocks, and they correspond to dark regions in the experi- 
mental schlieren. The light regions in both the constructed and experimental schlieren are identified 
with shear layers, and they are caused by density gradients near the surface. These shear layers lift 
off the surface when separation occurs, and the extent of the computed and experimental separation 
regions can be compared. The sharpness of the features in the constructed schlieren can be further 
improved by using grid adaptation to reduce solution errors. 

In figure 4, constructed and experimental shadowgraphs for a generic. National Aerospace Space 
Plane (NASP) nozzle configuration, the Single Expansion Ramp Nozzle (SERN), are shown. Both 
the freestream and nozzle flows are supersonic. The flowfield solution in figure 4(a) is a 2-D, ideal- 
gas simulation at the symmetry plane; the solution in figure 4(b) is 3-D (ref. 6). The interaction of 
the external flow with the nozzle plume is complicated, and all the main features observed in the 
experimental shadowgraph are captured in the constructed shadowgraphs. These features include the 
leading edge shock, separation from the nozzle, and the shocks, shear layers, and expansion fans 
caused by the interaction of the nozzle plume with the freestream. As with the schlieren, dark and 
light regions in the constructed shadowgraph correspond to dark and light regions in the experi- 
mental shadowgraph; however, in the shadowgraphs, each shock is described by the combination of 
a dark and light line. In the two constructed shadowgraphs, the effects of grid resolution and solution 
errors can be seen. The grids in the two 2-D solutions are very fine, and the constructed shadow- 
graph closely resembles the experimental shadowgraph. For the 3-D solution, although multiple 
grids were used, the gridding was much coarser than that used for the 2-D solutions. The shocks, 
shear layers, and expansion fans were not as well resolved in the 3-D flowfield solution and, hence, 
are much more diffuse in the shadowgraph. 



Constructed schlieren Experimental schlieren 


Figure 3. Constructed and experimental schlieren for a two-dimensional compression comer. Flow 
solution: Tokarcik et al. (ref. 4); experiment: Delery and Coet (ref. 5). 
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a) 2-D solution b) 3-D solution c) Experiment 

Figure 4. Constructed and experimental shadowgraphs for SERN. Two-dimensional flow solution 
and experiment: Ruffin et al. (ref. 6). 


The solutions in figures 3 and 4 are ideal-gas computations, and the refractive index is simply 
given by a constant times the density. For reacting flows, the refractive index is a function of the gas 
species. In figure 5, a constructed shadowgraph is shown for an axially-symmetric, real-gas, flow- 
field solution for a ballistic range shot of a hemisphere cylinder in a combustible mixture (refs. 7 
and 8). Also shown is an experimental shadowgraph (ref. 9). This flow is very complex, and many of 
the features observed in the shadowgraph are three dimensional and difficult to identify in contour 
plots. In the constructed shadowgraph most if not all of the features in the experimental shadowgraph 
are captured, and the computations and experiment can be compared on a one-to-one basis. The 
ability to construct shadowgraphs facilitates the comparison of computed and experimental results, 
especially for very complicated flows. 

In figure 6, a constructed and experimental infinite-fringe interferogram is shown for a ballistic 
range test of the Aeroassist Flight Experiment (AFE). The flowfield simulations were performed 
using an ideal-gas CFD solver (ref. 10); the grid was adapted to the flowfield solution using SAGE 
(ref. 1 1). A one-to-one comparison of the two interferograms is nof possible; the flowfield solution is 
for a model with no yaw, the experimental model did yaw. However, there are similarities. The con- 
structed and experimental interferograms both show similar light and dark patterns, and the number 
of fringe shifts are comparable. The location of the fringes provides information concerning the 
computed and experimental integrated index of refraction. 




Figure 5. Constructed and experimental shadowgraphs for a hemisphere-cylinder in a combustible 
mixture. Flowfield solution: G. Wilson (refs. 7 and 8); experiment: Lehr (ref. 9). 
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Constructed interferogram Experimental interferogram 

Figure 6. Constructed and experimental infinite-fringe interferograms for an AFE model. Three- 
dimensional flowfield solution: E. Venkatapathy et al. (ref. 10); experiment performed at NASA 
Ames Research Center’s ballistic range. 


CONCLUDING REMARKS 

The construction of interferograms, schlieren, and shadowgraphs from 2-D and 3-D flowfield 
solutions permits one-to-one comparisons of CFD and experimental results. Experimental interfero- 
grams provide quantitative information for integrals of the refractive index and density. This infor- 
mation can be used to verify 3-D flowfield solutions only when line of sight integrals are computed 
and interferograms are constructed from the flowfield solutions. The positions of shocks and shear 
layers are easily obtained from experimental schlieren and shadowgraphs. When contour plots are 
used to locate shocks and shear layers, interpretation of the contour lines is required. Furthermore, 
the 3-D nature of the flow is not shown. In constructed schlieren and shadowgraphs, the shocks and 
shear layers are defined by the same intensity variations as in the experimental images. Hence, com- 
parison of the positions of these flow features can easily be made. In addition, 3-D effects are 
included in the construction. 

CISS has proven its capability in constructing interferograms, schlieren, and shadowgraphs from 
a variety of ideal- and real-gas, 2-D, axially symmetric, and 3-D flowfield solutions. For the exam- 
ples shown here, the constructed and experimental images are very similar. However, the quality of 
the flowfield solution does affect the quality of the constructed images. Insufficient or improperly 
placed grid points can cause spreading of shocks and shear layers in the constructed images. CISS 
not only provides an excellent tool for comparing flowfield solutions with experiment, it can be used 
to identify problem areas in the flowfield solutions. 
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